Nonequilibrium critical dynamics of the three-dimensional gauge glass 
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We study the non-equilibrium aging behavior of the gauge glass model in three dimensions at the 
critical temperature. We perform Monte Carlo simulations with a Metropolis update, and correla- 
tion and response functions are calculated for different waiting times. We obtain a multiplicative 
aging scaling of the correlation and response functions, calculating the aging exponent b and the 
nonequilibrium autocorrelation decay exponent X c /z c . We also analyze the fluctuation-dissipation 
relationship at the critical temperature, obtaining the critical fluctuation-dissipation ratio Xoo. By 
comparing our results with the aging scaling reported previously for a model of interacting flux lines 
in the vortex glass regime, we found that the exponents for both models are very different. 

PACS numbers: 74.25.Qt,75.50.Lk,74.40.+k,05.70.Ln 



I. INTRODUCTION 

The possibility of a vortex glass phase in high-T c su- 
perconductors was proposed in 1989 by M. Fisher and 
coworkers J&L Throughout the years, two possible sce- 
narios have been discussed: either there is a continu- 
ous phase transition from the vortex liquid to the vortex 
glass at a finite critical temperature?^' 5 -' 6 -! 7 - or there is a 
crossover temperature 8 -'^'^ - below which the vortex liquid 
freezes in a glassy regime. While early theoretical^ and 
experimental studies 6 - supported the existence of a finite 
temperature critical point, more recent theoretical argu- 
ments tend to favor the freezing scenario for the vortex 
liquid at a crossover temperature. 8,9 Experimentally, it 
is difficult to differentiate between a critical point and a 
strong crossover in near equilibrium measurements in su- 
perconductors with random pinning JaIS The aim of this 
work is to show that the study of the nonequilibrium dy- 
namics can help to clearly distinguish among these two 
scenarios. To this end, we will study nonequilibrium cor- 
relation and response functions in a model for the vor- 
tex glass transition that has a finite temperature critical 
point (the unscreened gauge glass modelp and compare 
with previous results in a model that has a crossover 
freezing transition to a glassy regime (the interacting flux 
lines model) 

The gauge glass (GG) model^ 5 -*^ is a paradigmatic 
model for the vortex glass phase transitions^ It consists 
on the XY model with a random gauge potential vec- 
tor. The GG is assumed to represent the physics of a 
superconductor with disorder at relatively high magnetic 
fieldsS^ It has been found that this model has a finite 
T c in three dimensions 5 - ' 1 2 ' 13 and several numerical works 
have studied its critical behavior, considering both equi- 
librium properties as well as dynamical and transport 



properties 
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Early experiments* on 



the scaling of current- voltage curves in high-T c supercon- 
ductors obtained critical exponents v and z in reasonable 
agreement with the exponents obtained in simulations in 
the GG. However, in 1995 Bokil and Young 8 - showed that 
T c = in three dimensions when magnetic screening in 
a finite length scale A is included in the GG (by adding 



fluctuations of the vector potential) £ ' 19 ' 20 ' 21 Later ex- 
periments also showed that the current-voltage curves do 
not scale as expected?^ while other experimental groups 
report good scaling with the GG exponents^ Since all su- 
perconductors have a finite screening length A (the Lon- 
don penetration depth), the findings of Ref. 8 discard 
the scenario of a finite temperature critical point. Tak- 
ing into account this result, Zimanyi and coworker s 9 ' 22 
considered a model of interacting flux lines (IFL) with fi- 
nite interaction length A and in the presence of quenched 
disorder.— They found in this model that upon decreas- 
ing temperature the vortex-glass-like criticality is ar- 
rested at a crossover temperature signaling a freezing of 
flux lines. Below this temperature the time scales grow 
very quickly with a dynamical behavior similar to what 
is found in structural ("window") glasses^! 

This later result brings into attention the importance 
of the nonequilibrium behavior of vortex glasses as well 
as its possible relationship with the dynamical behav- 
ior of other glassy systems. In the recent years, there 
has been some progress in the general understanding 
of glasses through the study of their out-of-equilibrium 
dynamics^ 5 . A characteristic of relaxing glassy systems 
is the loss of stationarity reflected by their aging proper- 
ties, meaning that the dynamics of the system depends 
on the time elapsed after the preparation of the sample, 
t w £^ As a consequence, dynamic correlation functions 
dependend on two times, the "waiting time", t w , and the 
time t elapsed during the measurement. Also the lin- 
ear response functions show aging effects, being depen- 
dent on t w and t, and they are anomalous in the sense 
that they are not related to their associated correlation 
functions by the equilibrium fluctuation-dissipation theo- 
remi 25 i 26 The standard protocol for the study of aging in 
glassy system s 25 ' 26 consists on preparing the sample at a 
high temperature (T sta rt — ► oo) and then to quench it to 
a final low temperature Tf. At Tf, two-time correlation 
functions C(t, t w ) and response functions R(t, t w ) are an- 
alyzed, where the measurement starts for t >t w . In sev- 
eral systems a simple aging law with C(t, t w ) ~ C ag (t/t w ) 
is typically founds On the other hand, a "multiplica- 
tive" aging law of the form C(t,t w ) ~ t~ a C ag (t/t w ) has 
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been found in polymers in random media 27 and in sys- 
tems that are quenched at the critical temperature, like 
ferromagnetic spin model o 28 ^ 29 ' 30 ! 31 ' 32 ^ 33 and the Ising sp- 
ing glass i 34 i 35 

In the case of vortices in superconductors with 
quenched disorder, Bustingorry et alw^ have studied the 
nonequilibrium aging dynamics of the IFL model. They 
have found that it can be characterized by "multiplica- 
tive aging" , with a scaling form similar to what is found 
in polymers in random media.— In order to compare with 
this result, we will study here the the aging and nonequi- 
librium dynamics of the GG model at the critical tem- 
perature, calculating two-times correlation and response 
functions. Being at the critical point, a multiplicative 
aging law is expected, which in form should be similar 
to the type of aging observed numerically in the IFL. 24 
Our aim is to compare the aging exponents of the GG 
at T c with the ones obtained for the IFL in Ref. [24J. As 
we will show here, they turn out to be very different, 
a result that indicates that nonequilibrium aging exper- 
imental measurements could clearly distinguish among 
the two scenarios discussed in the first paragraph. 

The paper is organized as follows. In Section II 
we present the model Hamiltonian and the simulation 
method. In Section III we define the observables to be 
calculated. In Section IV we present our results for the 
correlation and response functions after a critical quench 
in the GG and we analyze their scaling with t w . Finally 
in Section V we discuss our results comparing with the 
disorder-free 3D XY model, the 3D ISG and the IFL. 



are analyzed, which depend on both the waiting time t w , 
when the measurement begins, and a given time t > t w . 

For the Monte Carlo simulation local changes in the 
phases Oi — * 9[ are accepted with probability given by 
the Metropolis rate 



p(0i -» e[) = min{l,exp(-/3A#)}. 



(2) 



Here /3 is the inverse temperature and AH is the energy 
difference corresponding to the proposed phase change. 
In equilibrium simulations, the acceptance window for d' i 
is usually chosen less than 2n and dependent on tem- 
perature in order to optimize the updating procedure. 
Because we are interested in studying a nonequilibrium 
process, for local phase changes we will use the full 2ir 
acceptance angle window for the new phases 0'., follow- 
ing the criterion used in Ref. [l7lfT8l . This is done in order 
to avoid the possibility that a limited acceptance angle 
might introduce an artificial temperature dependence in 
the relaxation. 



III. 



OBSERVABLES 



For simplicity, we will focus on the study of the two- 
times autocorrelation function defined as 



C(t, t w ) 



N 



N 



5^cos[0<(t) -Oi(t w )] 



(3) 



II. MODEL AND MONTE CARLO 
SIMULATIONS 



The hamiltonian of the three-dimensional 
model- is given by 



(3D) GG A 



H = 



-J cos( 



6 4 - A 



1,3) 



(1) 



where Oi represents the superconducting phase at site i, 
and we sum over nearest neighbors on a cubic lat- 
tice of linear size L (N = L 3 ). The A^.j are quenched 
random variables uniformly distributed in the [0, 2ir] in- 
terval (Aij = —Aj t i); and J is the coupling between 
nearest neighbors. The phases Oi can be represented as 
the angle of classical two-dimensional spins of unit length, 
Si = (cos 9i, sxnOi). In this work periodic boundary con- 
ditions are applied in the phases 6i. The energy is nor- 
malized in units of J, temperature in units of J/fcs, and 
time scales are measured in number of Monte Carlo steps 
(for full sweep for the iV sites). 

The out-of-equilibrium protocol used in this work, con- 
sists on a quench at time t = from a T = 00 state to 
the critical temperature T c . For the 3D GG model, the 
critical temperature is T c = 0.46(1), according to Ref. [l7l. 
From this initial condition different two-times quantities 



where (...) indicates an average over different thermal 
histories (different initial configurations and realizations 
of the thermal noise), and [...] a v represents a disorder 
average over different samples (different realizations of 



h3 



The corresponding two-times linear autoresponse func- 
tion is 



R(t, t w ) 



N 



I N 

E 



SSi(t) 



— ^ 5hi(t w ) 



(4) 



In this case, the system is perturbed by applying an 
infinitesimal external field hi conjugated to Si. This 
corresponds to a perturbing term in the hamiltonian 
AH = — hi - Si. In numerical simulations it is more 
convenient to calculate the integrated responses: either 
the thermoremanent response MxRM(t,t w ) or the zero- 
field-cooled response MzFc{t, tw), which are obtained by 
switching on the perturbation only for times t <t w and 
t > t w , respectively. If we define the reduced integrated 
responses by 



T 

PTRM(t,t w ) = —MTRM(t,t w ) 
h 

T 

PZFc{t,tw) = -rMzFc(t,t w ), 



(5) 
(6) 
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then we can relate these functions to R(t,t w )£2- 



Ptrm (t, t w )=T / duR(t, u) 



PzFc(t,t w ) = T I duR(t,u). 



(7) 
(8) 



If the perturbing field in each site, hi, is random and its 
two components are independently drawn from a bimodal 
distribution ±h, both reduced integrated responses are 
given by 3 -^ 



p(t, t w ) 



T 



N 



h i ■ S i 



(9) 



In this work we show simulation results for systems of 
size L = 30 at the temperature T c = 0.46. We take the 
disorder average over 60 samples and, for each sample, 
we carry out a number of 10 thermal histories. The mag- 
nitude of the external fields used for the calculation of 
the response function was h = 0.1 and h = 0.05. 



IV. CRITICAL QUENCH 

At the critical temperature T c the equilibrium auto- 
correlation relaxation time r e increases with the size as 
t c ~ L Zc , where z c is the dynamical critical exponent. 
For simple ferromagnetic systems, it is known that in a 
critical quench spatial correlations over a length scale of 
I ~ t 1 / Zc are just as in the equilibrium critical state. This 
means that the system appears critical on scales smaller 
than I, while it appears disordered on larger scales.— 

Considering this, the autocorrelation function © is 
expected to behave a a 29 i 35 



C(t, t w ) =t w fc[ — 



(10) 



For r <C t w (t = t — t w ), the critical scaling function 
fc{x) behaves as 

f c {x)~[x~l]- h . (11) 
For d dimensional spin glasses it has been foun d 34 ! 36 that 

d- 2 + 7/ 



b = 



2Zr 



(12) 



Here r] is the static critical exponent associated to the 
pair correlation function. Then, for r <C t w the autocor- 
relation function is 



C(t,t w )~T- b . (13) 
On the other hand, when r 3> t w we have: 

fc(x)~A c x- Xc ' z ^ (14) 
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FIG. 1: (Color online) (a) Autocorrelation function Cvst for 
eight different waiting times as indicated, (b) Data collapsing. 
Inset: autocorrelation function C vs t/t w . Time scales are 
measured in number of Monte Carlo steps (for full sweep for 
the N sites). 



where Ac is a constant and Ac is known as the autocor- 
relation exponent^ 4 - 

The Fig.QJa) shows the behavior of the autocorrelation 
function C(t,t w ) vs r for eight different waiting times 
t w . For t <§C tw we can fit a power law behavior, t~ 6 , 
with exponent b w 0.06, as shown in the plot. For r > 
t w we observe that the correlation function shows aging 
since it depends strongly on the waiting time t w . In the 
inset of Fig.Q] (b) we show the result of assuming "simple 
aging" by plotting C as a function of the ratio t/t w . We 
observe that the data does not show a good collapse in 
a single curve in this case. A good data collapse can 
be obtained if the "multiplicative aging" scaling form of 
Eq. lfTO")) is assumed. In Fig. Q] (b), we show the plot of 
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t b w C vs. t/t w , obtaining good scaling with b = 0.06(1). 
This value is in good agreement with Eq. (JT2J) : if we take 
from previous simulations of the 3D GG model the values 
T) = -0.47(2)" and z c = 4.5(1)±£ we obtain b = 0.059(4). 
We now consider the case of r 3> t w . In the Fig. Q] (b) 
we show that in this case the asymptotic behavior of C 
has a power law form consistent with Eq. (fT4"]) . We obtain 
A c = 1.55(5) and \ c /z c = 0.49(2). 

We now consider the autoresponse function. In this 
case the scaling equation is^ 

R(t,t w )=t w - i - a f R (j-y (15) 

For r ^> t w is expected that 

f R (x) ~ A R x^^ 

where Ar is a constant amplitude. Then, if we con- 
sider Eq. l[7]). we obtain the thermoremanent reduced in- 
tegrated response 

pTRM(t, t w ) = t w - a f p (j-^j , (16) 

where 

f p (x) ~ A p x- X *l*°, (17) 

and A p is a constant. For critical systems it is expected 
that a = b, and Ac = A#. 

In Fig. [5] (a) we show the behavior of Ptrm function vs 
r for different waiting times and h = 0.05, observing that 
the response function also shows aging. In the inset of 
Fig. [2](b) we show Ptrm vst/t w : no data collapse is ob- 
served. In Fig. [U(b) we plot V^ptrm vs t/t w , obtaining 
good data collapse for a = 0.06(1). The asymptotic be- 
havior observed for large temporal separations (r 3> t w ) 
can be fitted with a power law with A p — 0.18(2) and 
\ R /z c = 0.52(2). 

At thermodynamic equilibrium the fluctuation- 
dissipation theorem (FDT) is satisfied, and therefore 
it is possible to draw a simple expression to relate the 
correlation with the thermoremanent reduced integrated 
response 

PTRM{t-t w ) = C{t-t w ). (18) 

Then, for a parametric plot of ptrm vs. C, we should 
obtain a straight line of slope +1 in equilibrium. On the 
other hand, for a nonequilibrium process the FDT is not 
fulfilled, and it has been proposed a generalized relation 
of the form 

pTRM{t,t w ) = X(t,t w )C(t,t w ), (19) 

where X(t,t w ) is called the fluctuation-dissipation ratio 
(FDR)i2I For 1 <g; t w -C t we can estimate the limit : 

Xco = lim lim X(t,t w ). 

t w — >oc t — >oo 
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FIG. 2: (Color online) (a) Function ptrm vs t for different 
waiting times and h — 0.05. All symbols are the same as in 
Fig.[T](a). (b) Data collapsing. Inset: Function ptrm vs t/t w 
for different waiting times and h — 0.05. 

For a critical quench, Xaa is expected to be universal in 
the sense that it does not depend neither on the initial 
conditions nor on the details of the dynamics.— 

The Fig.[3]shows the t-parametric plot of PTRAt(t,t w ) 
vs C(t,t w ) for different waiting times t w . In the quasi- 
cquilibrium regime (r <C t w ) the FDT holds and X = 1. 
This is observed in Fig. [3] where for large t w the slope 
in the FDT plot tends to 1. On the other hand, when 
1 <C t w <C t, we observe that X < 1 and FDT is violated. 
From Eqs. (JTOj), (O, (JTSJ), ([TT |) and ([19 ]) . it is possible 
to demonstrate that 

Xoo = A p /A c . 

From the amplitudes Ac and A p calculated previously 
we obtain X^ = 0.12(2). 
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FIG. 3: (Color online) FDT plot for different waiting times 
and h = 0.05. All symbols are the same as in Fig. [T](a). 



We can obtain Xoo numerically from the FDT plot of 
Figure 3 in two different ways. First, in Fig.|4](a) we show 
a data collapse plotting ptrm (t t t w )t%, vs C(t, t w )t b w . The 
slope near the origin corresponds to the limit of FDR, 
obtaining X^ — 0.12(1). Second, in Fig. 0] (b) we plot 
the ratio t w a pTRM/t w b C vs t w b C for all times t and t w . 
We see that the ordinate at the origin is compatible with 
Xoo ~ 0.12. 



V. DISCUSSION 

We have therefore obtained the exponents that charac- 
terize the aging scaling of the 3D GG model at the criti- 
cal temperature, obtaining b — 0.06(1), Xc/z c — 0.49(2), 
a = 0.06(1), \ R /z c = 0.52(2), and X^ = 0.12(1). A 
comparison with results in the corresponding model with- 
out disorder, i.e. the 3D XY model, shows that the 
exponents are very different, as expected. Abriet and 
Karevski^ have found a = b = 1/2, Xc/z = Xr/z 
I ..'!-">. which shows that in the disorder-free case the "mul- 
tiplicative" scaling of aging is stronger (i.e., much larger 
values of a and b). Also a very different limit of the 
FDR, = 0.43(4), was obtained, as expected since 
they should belong to different universality classes. 

On the other hand, the critical aging of the GG 
turns out to be very similar to results found in the 3D 
ISGi 34 i 35 It has been found for a binary distribution in 
the exchange couplings: a — 0.060(4), b = 0.056(3), 
X c /z = 0.362(5), Xr/z = 0.38(2) and X x = 0.12(1); 
and for a Gaussian distribution in the exchange cou- 
plings: a = 0.044(1), b = 0.043(1), X c / z = 0.320(5), 
Xr/z = 0.33(2) and = 0.09(l)^i As we can see, the 
aging exponents are slightly smaller in the 3D ISG (mean- 
ing a weaker multiplicative aging), but the values of X^ 
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FIG. 4: (Color online) (a) Data collapsing for all curves in 
Fig-E] (b) Plot of ratio t w a p T RM/t w b C vs t w b C. All symbols 
are the same as in Fig. [T] (a). 



are comparable within the statistical error with the GG. 
It has been argued in Ref. 28.29.30 that X^ is a novel 
universal quantity of non-equilibrium critical dynamics. 
It is interesting to observe that spin glass models like the 
ISG and GG have similar values of X^. However, let 
us remind that most work on the equilibrium properties 
of the GG and the ISG show evidence that they are in 
different universality classes . 5 ' 11 ! 12 ' 15 ! 18 

As discussed in the introduction, two types of behav- 
iors are predicted for the vortex glass transition: either 
a finite temperature critical point as described by the 
GG, or a crossover temperature as found in the IFL. 
The work of Bokil and Young 8 discards the first sce- 
nario. However, from the experimental point of view, 
the issue is not clearly resolved since several exper- 
iments show reasonably good scaling at a transition 
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temperature^ Comparing our results with the behavior 
found by Bustingorry et al~ in the IFL, we believe that 
a study of the non-equilibrium dynamics of the vortex 
glass transition will be able to clearly distinguish be- 
tween the two models. In simulations of the IFL , off- 
equilibrium correlation and response functions arc found 
to scale as B(t,t w ) ~ t™B(t/t w ), and also it is found 
that B{x) ~ x a for x > 1^ The first thing to no- 
tice when comparing with critical aging in the GG is 
that, in the IFL, the corresponding multiplicative aging 
exponent b and the autocorrelation exponent Xc/ Z co- 
incide b = Ac I z — a (and similarly for the response 
exponents, a — Xr/z = a). The second main differ- 
ence is that the "multiplicative" aging in the IFL is very 
strong, &ifl = ol ~ 0.25 — 0.35 when compared with the 
GG at the critical point, 6qg ~ 0.06. In the IFL the 
competition between the elasticity of each flux line and 
the randomness of the pinning potential leads to a very 
slow dynamics and strong aging with a large a exponent. 
Moreover, a depends with the disorder strength in the 
IFL, being larger for stronger disorder. 

In this work we have analyzed the correlation and in- 
tegrated response functions defined in Eqs.© and ©, 
since they are easy to compute numerically. In exper- 
iments in superconductors, correlation (response) func- 
tions of other observables like magnetization (magnetic 



susceptibility) and voltage (resistance) are, of course, 
more available. Aging scaling similar to Eqs. tflT))) and 
(fTo) will also be expected in these other physical quanti- 
ties. For example, one possible way of observing exper- 
imentally the out-of-equilibrium dynamics of the vortex 
glass is monitoring the magnetic relaxation as done in 
Refl38l for granular samples. Another possibility is to per- 
form electrical transport experiments near the transition 
temperature following a protocol similar to the one used 
by Ovadyahu and coworkers to study the electron glass 
regime in Anderson insulators^ In either case, an analy- 
sis of the multiplicative aging exponent will clearly distin- 
guish among the two different models (6ifl ~ 0.25 — 0.35 
or b GG » 0.06). 
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